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Abstract 

The behaviour of a model glass forming liquid is analyzed for a range of 
densities, with a focus on the temperature interval where the liquid begins 
to display non-Arrhenius temperature dependence of relaxation times. Ana- 
lyzing the dynamics along with properties of local potential energy minima 
sampled by the liquid, a crossover or onset temperature T s is identified below 
which the liquid manifests slow dynamics, and a change in the character of 
typical local potential energy minima. 
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The phenomenology of glass formation, whereby a liquid cooled to low temperatures 
manifests increasingly sluggish dynamics, and at low enough temperature transforms into 
an amorphous solid, has come under increasing study in recent years [0,0. In addition to 
the glass transition temperature T g where the liquid transforms to a glass in the laboratory, 
much attention has been devoted to the Kauzmann or ideal glass transition temperature Tk 
||-|9|], which locates a putative thermodynamic phase transition underlying the observed glass 
transition, and the mode coupling critical temperature T c [p^ -|l2|l, which is now commonly 
viewed as marking the crossover to a regime of 'activated dynamics' discussed long ago by 
Goldstein ||13|| . Relatively less attention has been paid to the range of temperatures (or 
equivalently, a crossover temperature) where the nature of the dynamics changes from what 
may be termed as 'normal' to what is commonly referred to as 'slow dynamics' fl4HT6||. The 



purpose of this letter is to demonstrate, on the basis of computer simultions of a model 
liquid, that such a crossover temperature can be identified by various indicators that mark 
a change in the nature of local potential energy minima the liquid samples, in addition to 
changes in the dynamics itself. 

Results presented here are from molecular dynamics simulations of a model liquid which 
is a binary mixture consisting of 204 type A and 52 type B particles. The particles interact 
via the Lennard- Jones (LJ) potential, with parameters €ab / £aa = 1-5, cbb/caa = 0.5, 
ctab/vaa = 0.8, and o~bb/ ® aa = 0.88, and m B /mA = 1. This system has been extensively 
studied as a model glass former fLT , pL7|, 14| ,|6|,p|Jl^| . The LJ potential is modified with a 
quadratic cutoff and shifting at r^f = 2.5<j a p [[HJ. All quantities are reported in reduced 
units, length in units of uaa-, volume V in units of a\ A (density p = N/V, where N is 
the number of particles, in units of a AA = po), temperature in units of £aa/ energy 
in units of eaa and time in units of r m = (a 2 AA m/ 6aa) 1 > where m = tua = wig is the 
mass of the particles. Molecular dynamics simulations are performed over a wide range 
of temperatures at each density, and local energy minimizations are performed for 1000 
(ksT/tAA < 1.) or 100 {ksT jtAA > 1-) configurations to obtain typical local energy minima 
or 'inherent structures' [^] sampled by the liquid. Inherent structures are obtained for a 



subset of temperatures at intervals that vary between 0.2 to 0.4 near the crossover T, which 
thus limit the precision of estimation, in addition to the inherent lack of sharp definition of 



a crossover temperature. Further details concerning the simulations may be found in |18 
where a study of the dynamics and thermodynamics of this model was reported over a wide 
range of densities and temperatures. 

Data from these simulations are analyzed here for five densities to study the crossover 
temperature where the liquid begins to manifest slow dynamics. The dynamics is probed 
by calculating the self intermediate scattering function F s (k,t) for the A particles, which is 
the Fourier transform at wavevector k of the van Hove self correlation function 

i N A 

Gs(r,t) = —Y,(K\r l (t)-r(0)\-r)), (1) 

"A i= i 

which describes the probability of finding at time t a particle at distance r from its location 
at t — 0. At each density, the wavevector studied is close to the first peak of the structure 
factor. Relaxation times r(T) are obtained as a function of temperature T, at each density, 
from F s (k, t) with the definition, F s (k, r) = 1/e. At high T, r display Arrhenius behavior (i. 
e., t = T exp(E /kBT) where Tq and E are constants) as shown in Fig. 1(a). Deviations 
from Arrhenius behavior can be seen clearly by defining a T-dependent 'activation energy' 

E(T) = k B T\og(r/r ) (2) 

which has the constant value E for T where the Arrhenius form is valid. High T values of r 
are used to estimate Tq and E which are then used to transform the full data set, following 
the procedure in |14[]. The resulting E(T) data (normalized to Eq) are shown in in Fig. 



1(b), which display for all densities p the crossover to non- Arrhenius behavior, similar to 
the finding in |TJj] for p = 1.2p . The crossover T identified this way are shown in Fig. 2. 
Also shown are the ideal glass transition line ('Tr-'; from [18[]) and the VFT temperature T 
obtained from diffusivity data in ||18|| . The mode coupling temperatures T c have not been 
estimated as a function of density in this study, but will be presented in a future article. 
At p = 1.2p , the value of T c estimated in |TTJ is ksTc/e^A — 0.435, while from the present 



study, kBT s /eAA = 1-0 and kBTx/eAA = 0.2976. Therefore, the mode coupling T c at this 
density is less than half the onset temperature T s , and higher than the ideal glass transition 
temperature Tk by roughly 50%. 

Properties of potential energy minima sampled by the liquid show a corresponding 
crossover behavior, which is now discussed. A number of recent studies [pi],P-P| J22] , |25| j9] jr8|j2^] 
have focussed on analyzing the thermodynamics of glass forming liquids on the basis of a 
decomposition of the configurational space of the liquid into basins, as a route to estimating 
the configurational entropy, which is given by the multiplicity of such basins. In the inherent 
structure approach [pOfl , basins are identified as basins of the local potential energy minima 
(inherent structures), which permits the partition function to be written as (details may be 



where <fr is the potential energy per particle of the inherent structure, Q((j))d(j) is the number 
of inherent structures in the range (0, 4> + d(p), exp(— j3 fbasini^i T)) is given by the restricted 
partition function integral over the basin of an inherent structure of energy 0. The configu- 
rational entropy density is defined as <S C (0) = % lnf2(0) and is related to the T dependent 
configurational entropy by S C (T) = f d<p S c (<p)P(<f),T). The probability of sampling an 
inherent structure of energy <ft is given by 



The distribution of inherent structures sampled is calculated directly from simulations, by 
performing local energy minimizations for a sample of equilibrated configurations. In addi- 
tion, the total free energy of the system (and hence Qn) is calculated from thermodynamic 
integration as outlined in [0. Knowledge of the basin free energy f basin is then sufficient 
to invert the above relationship to obtain the configurational entropy density S c . One must 
obtain the same estimates of S c from P(4>, T) at different temperatures, in the range of 
values where the corresponding P(<p, T) overlap. In previous work ||[7],[24|] it has been 



found in @g0) 




(3) 



P(0, T) = exp (-/3N(4> + f basin - S c )) /Q N {p, T). 



(4) 
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observed that one indeed obtains overlapping S c estimates from different temperatures at 
sufficiently low temperatures. Results in these analyses was based on evaluating the basin 
free energy with a harmonic approximation, wherein the basin free energy is given in terms 
of the frequencies V{ which are related to the eigenvalues A, of the Hessian matrix (matrix 
of second derivatives of the potential, defined at the energy minima) by 27rz/j = The 
basin free energy is given in this approximation by: 

1 3N hvi 

Pf basin = J-f J2 111 ( 5 ) 

The eigenvalues Aj are obtained numerically for the sample of inherent structures obtained at 
each temperature and density. Estimates of S c from different temperatures for each density 
are shown in Fig. 3. It is seen that at each density, the S c estimates from low temperatures 
overlap with each other, but estimates from higher temperatures fail to do so (each T is not, 
however, individually labeled in the figure). The harmonic approximation for describing the 
properties of basins therefore fails at high temperatures, indicating a qualitative change in 
the nature of the basins sampled at temperatures below and above a crossover temperature. 
The crossover temperatures estimated in this fashion are shown in Fig. 2. The criterion 
used here, as well as in Fig. 4, 5, is stated in the figure caption. 

It has been observed |7|, ^5| , ^| that distribution of inherent structures fl(4>) is well de- 
scribed by a Gaussian (equivalently, the configurational entropy density S c by an inverted 
parabola). Figure 3 shows fits to S c at each density of the form, 

SM)lhB = a - (6) 

The parameters a, 0o and a depend on the density p but not on the temperature T. Based 
on the assumption of a Gaussian fl(4>) and a constant (with respect to <fi) Pf basin in Eq. 



(0) it has been observed p5|,|7| that the average inherent structure energy <j>{T) should vary 
with temperature as <f>(T) = O — o 2 jlksT . However, the dependent component of fifbasin 
(arising from the dependence of the frequencies z/, on the basin occupied, indexed by its 
energy) is not constant but varies roughly linearly with ]27J (see also P5| , PP| ). Hence, the 
form 



/3/ 6as «(0) = /3/o + W-0 O ), (7) 

where f3fo has no dependence, is employed to describe the T and (j) dependence of (5 j 'basin- 
Inclusion of this <fi dependence does not alter the prediction of the 1/T dependence but 
redefines the quantitative expression to 

0(T) = - a 2 /2k B T, (8) 

where the effective mean IS energy = 0o — ^fo~ 2 /2. It has been shown elsewhere [27 



that Eq. (||) indeed describes well the average inherent structure energies, at low temper- 
atures, calculated directly from simulations. However, at high temperatures, one observes 
deviations from the 1/T behavior (see also f24j|). A scaled plot of 4>{T) — </>q * vs. a 2 /2ksT 
is shown in Fig. 4. In order to identify the crossover temperature as accurately as possible, 
the fit parameters involved, a, 0o> o~ and 5f have been estimated by fitting jointly the low 
temperature data sets of S c (4>) (Eq.(|6])), fifbasin (Eq.(^)) and 0(T) (Eq.(|])). The fit pa- 
rameters are shown in Table I. Crossover temperatures where <f)(T) exhibits deviations from 
the 1/T form are shown in Fig. 2. The analysis here makes more precise, and quantita- 
tive, the observation in |TJj of a correlation between the onset of non-Arrhenius relaxation 
and a change from nearly constant to T dependent values of the average inherent structure 
energies. 

When inherent structure basins sampled by the liquid are well described by a harmonic 
approximation, one must expect that the potential energy per particle u{T) of the equili- 
brated liquid is related to the inherent structure energy 0(T) by u{T) = 4>{T) + 3/2/c B T. 
Thus significant differences between the quantity Au(T) = u(T) — 3/2fc#T and 0(T) also 
indicate anharmonicity of the basins sampled |30,24j. Fig. 5 shows Au(T) plotted against 



0(T) for each of the studied densities. The range of inherent structure energies sampled is 
non-monotonic in density, with the lowest energies obtained at p = 1.2p . For each density, 
there is a low energy range for which the harmonic expectation (indicated by the solid line) 
is met to a good extent (the deviations from the harmonic expectation in this range is less 
than 4% in all cases). At higher energies, it is seen that a fairly sharply defined energy (or 



temperature) exists where the two quantities deviate significantly. Such deviation permits 
the identification of a crossover temperature. The crossover temperatures so obtained are 
plotted in Fig. 2. 

Comparison of the various crossover temperatures in Fig. 2 identified from characterizing 
the dynamics (Fig. 1) and the inherent structures sampled (Fig.s 3, 4, 5) show that the 
crossover from 'normal' to 'slow' dynamics is associated with a change in the nature of 
local potential energy minima the liquid samples. The crossover is particularly broad for 
p = 1.35 which correspondingly shows the poorest agreement between different estimates. 



The correspondence of changes in IS properties with liquid structure was studied in [16 



and more recently in [29| for water. Some analysis of this crossover has been made for 



spinglass models, e. g. [^,[3T]]. In |23| it is shown that finite range spin glass models show 



behavior of average energies of local minima which are similar to liquids. Similar findings 
for the nearest neighbor ±J spinglass model have been reported in |31|. Analysis of the 
free energy landscape for the hard sphere fluid |32[ identify a crossover density where glassy 
free energy minima come into existence, which may be related to the crossover discussed 
here. Frustration limited domain theory |I3J discusses an avoided critical point temperature, 
which has been identified with the crossover temperature discussed here. 

Nevertheless, the precise nature of the crossover is at present unclear. Although Fig.s 3-5 
indicate a change from 'harmonic' basins at low temperature to those that are anharmonic at 
T above the crossover T, detailed analyses 0,^,0 indicate that at temperatures between 
mode coupling T c and the crossover temperature T s , the dynamics appears via motion along 
directions in configuration space with negligible energy barriers, the barriers being instead 
entropic ||32|| . Therefore, even though the dynamics in the range of T c to T s is 'landscape 
influenced' it is 'landscape dominated' only below T c Notwithstanding the quanti- 
tative measures presented here, the precise sense in which the harmonicity of local minima 
changes across T s is unclear. Very little theoretical effort has been directed at understanding 
the nature of the crossover at T s . Results presented here, along with related evidence cited, 
call for a focus on phenomena pertaining to the onset temperature of slow dynamics, in 



addition to more thoroughly studied phenomena at lower temperatures. 



p/Po 
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4>o/ € AA 






1.10 


1.020 


-6.579 


0.312 


0.616 


1.15 


0.963 


-6.680 


0.379 


0.603 


1.20 


0.921 


-6.700 


0.470 


0.455 


1.25 


0.875 


-6.642 


0.550 


0.350 


1.35 


0.860 


-6.080 


0.870 


0.159 



TABLE I: Fit parameters quantifying the distribution of inherent structures, and the depen- 
dence of vibrational free energy on the inherent structure energy, described in Eq.s (Eq.flf^)) 
and (Eq.©). 
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FIG. 1. (a) Arrhenius plot of relaxation times for different densities. At each density, a crossover 
from Arrhenius behavior at high temperature to super-Arrhenius behavior is observed, (b) Tem- 
perature dependent 'activation energies' E(T) normalized to the high temperature value Eq. 
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FIG. 2. Crossover temperatures estimated from (i) relaxation times(Fig. 1), (ii) Configurational 
entropy density S c (4>) (Fig. 3), (iii) T dependence of IS energies (Fig. 4) and (iv) the energy 
difference between instantaneous and IS configurations (Fig. 5). 
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FIG. 3. Configurational entropy density S c (4>) estimated as described in the text 
the (lower) value of <j> where S c (4>) = from Eq.(g). At each density (each overlapping set of 
curves), S c (cp) estimates obtained from a range of temperatures are shown (not individually la- 
beled). T* = ksT/eAA in the legends. It is seen that at each density, a number of data sets 
overlap with each other to a very good extent (these correspond to low temperatures) while some 
data sets show much poorer overlap (these correspond to high temperatures). The failure of high 
temperature S c {<p) curves to overlap with those at lower T signals a crossover in behavior. Data 
points for the lowest temperature where such lack of overlap is discernible are displayed using filled 
symbols and labeled. These temperatures provide estimates of the crossover temperature. 
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FIG. 4. Scaled plot of the T dependence of IS energies described in the text. Deviations from 
the expected linear behavior at high T marks the crossover. The crossover temperature is identified 
as the lowest T for which the IS energies show significant deviation from the low T linear behavior. 
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FIG. 5. The equilibrium potential energy u(T) minus 1.5ksT (the harmonic expectation) is 
shown against IS energies. When the harmonic approximation is valid one expects the data to fall 
on a straight line of unit slope. At each density the temperature where u{T) — l.bksT shows a 
maximum before dropping below the reference line is identified as the crossover temperature. 
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